function popm=Mutate(pop,lb,ub)
    nx=numel(pop);
    change=randperm(nx,1);%随机得到变异点位
    popm=pop;
    x1min=lb(1);x2min=lb(2);x3min=lb(3);
    x1max=ub(1);x2max=ub(2);x3max=ub(3);
    if change==1
        popm(1,change)=x1min+(x1max-x1min).*rand();
    elseif change==2
        popm(1,change)=x2min+(x2max-x2min).*rand();
    elseif change==3
        popm(1,change)=x3min+(x3max-x3min).*rand();
    end
end